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ABSTRACT 


A  real  time,  man  controlled  simulator  exhibiting  the 
characteristics  of  a  3000  Ton  Surface  Effect  Ship  (3K-SES) 
was  developed  using  a  set  of  simplified  5  degree  of  freedom 
equations  of  motion. 

A  software  program  was  developed  to  compute  and  present 
to  the  SES  pilot  on  a  real  time  graphics  output  "alpha- 
numeric-perspective" visual  cues  pertaining  to  the  vehicle 
states. 

Man  generated  control  responses  during  3K-SES  accidental 
effector  nozzle  failure  tests  and  the  requirement  for 
computer  assist  control  systems  during  the  recovery  phase  of 
such  failures  was  examined. 
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I.   INTRODUCTION 

Extensive  analysis  of  the  dynamics  of  the  Captured  Air 
Bubble  (CAB)  Surface  Effect  Ship  (SES)  has  been  accomplished 
in  the  past.   Such  analysis  has  produced  equations  of  motion 
of  the  SES  type  vehicle  in  six  degrees  of  freedom  [Ref .  1] . 
Analysis  at  the  Naval  Postgraduate  School  has  centered 
around  the  XR-3  and  100-B  SES  [Refs.  2,3],  with  current 
analytical  effort  being  directed  toward  the  3K-SES  [Ref.  4] . 
The  primary  tool  for  study  of  the  3K-SES  characteristics 
is  the  highly  comprehensive  five  degree  of  freedom  (DOF) 
SES  DBS  Program  developed  by  the  David  W.  Taylor  Naval  Ship 
Research  and  Development  Center  [Ref.  5] .   A  version  of 
this  program  used  at  NPS  is  called  DBSIM5D.   This  data  base 
program  allows  very  detailed  study  to  be  conducted  on  the 
3K-SES  response  to  precomputed  control  inputs  such  as  step 
or  ramp  thruster  commands  and  their  associated  accidental 
failure  possibilities. 

A  natural  extension  in  the  study  of  control  of  the 
3K-SES  is  to  evaluate  the  vehicle's  response  to  an  on-going 
"man  generated"  control  input  on  a  real  time  basis.   The 
rationale  for  developing  a  "man-in~the-loop"  SES  simulator 
would  be  to  1)  test  the  man  in  responding  to  accidental 
failure  situations  that  may  be  experienced  during  opera- 
tional conditions,  2)  train  personnel  in  operator  control 
techniques  unique  to  SES  type  vehicles,  3)  provide  a  means 
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to  test  hardware  to  be  utilized  in  SES  craft,  and  4)  provide 
a  design  tool  for  development  of  control  schemes  or  system 
structural  changes  to  improve  performance. 

Accordingly  a  simplified,  non-linear  five  degree  of 
freedom  Real  Time  Simulator  (RTS5D)  for  the  3K-SES  was 
developed  such  that  a  continuous,  on-going,  observable, 
real-time  solution  is  interfaced  with  man-generated  control 
inputs . 
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II.   EQUATIONS  OF  MOTION 

A.   COORDINATE  SYSTEMS  AND  ASSUMPTIONS 

It  was  recognized  that  the  RTS5D  would  require  accurate 
and  yet  simplified  equations  of  motion  when  compared  to 
those  of  the  more  extensive  DBSIM5D.   The  rationale  for 
requiring  simplification  was  to  start  with  a  right  hand 
side  (RHS)  that  could  be  utilized  effectively  in  a  real 
time  simulation.   To  satisfy  this  requirement  a  model  was 
selected  that  would  use  point  source  of  drag  forces  con- 
sisting of  only  a  few  lumped-coefficient  terms.   The 
equations  used  are  based  upon  those  developed  in  Ref.  1 
which  follow  this  concept.   Additional  assumptions  used  in 
the  simplified  5  DOF  equations  were  as  follows: 

1)  All  accelerations  are  measured  at  the  center  of 
gravity. 

2)  Vehicle  is  "f ree-to-heave" . 

3)  All  cross  coupled  moments  of  inertia  are  zero. 

4)  All  cross  products  of  angular  velocities  in  force 
equations  are  zero. 

5)  All  roll  and  pitch  angles  used  in  force  equations 
are  subject  to  small  angle  approximations. 

6)  Calm  water  conditions. 

7)  Mass  and  mass  distribution  of  vehicle  are  constant. 

8)  Effect  of  changes  of  aerodynamic  force  are  neglected, 
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utilizing  the  above  assumptions,  the  simplified  5  DOF 
equations  of  motion  are  developed  using  the  coordinate 
systems  shown  in  Fig.  1  and  Fig.  2,   the  force  vector 
diagrams  shown  in  Fig.  3,  Fig.  4A,  Fig.  4B,  and  the 
following  equations  from  Newton's  laws  of  motion: 


SURGE 

m(u  - 

vr) 

^ 

X 

SWAY 

m  (v  - 

ur) 

= 

Y 

YAW 

^z  ^ 

= 

N 

PITCH 

1^4 

= 

M 

ROLL 

^x^ 

= 

K 

where : 
m 
u 

V 

r 
P 

q 

i„ 


"X 


X 
Y 

N 


mass  of  the  rigid  ship 

velocity  in  x-direction  (SURGE) 

velocity  in  y-direction  (SWAY) 

angular  velocity  about  the  z-axis 

angular  velocity  about  the  x-axis 

angular  velocity  about  the  y-axis 

moment  of  inertia  about  the  z-axis 

moment  of  inertia  about  the  y-axis 

moment  of  inertia  about  the  x-axis 

summation  of  forces  in  x-direction 

summation  of  forces  in  y-direction 

summation  of  moments  about  the 
z-axis 


slugs 

ft/sec 

ft/sec 

deg/sec 

deg/sec 

deg/sec 

lb  -ft^ 
m 

lb  -ft^ 
m 

lb  -ft^ 
m 

lb, 

ft-lb. 
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M       =      suimnation  of  moments   about   the 

y-axis  ,      ft-lb^ 

K   =   summation  of  moments  about  the 

X-axis  ,   ft-lb^ 


Additionally  the  following  navigation  relationships 
are  utilized  (see  Fig.  1) , 


X  =   u  cos  i|i  -  V  sin  ^ 


Y  =   u  sin  \b   +  V   cos  \b 
o 


2  2 

V  =    u   +  V 

s 


S  =   tan'""'"  (-v/u) 


where:       4;   =  heading  angle 
3   =  drift  angle 
6   =  thrust  vector  angle 

Also,  X        =/xdt  u      =      /udt 

o  •'      o  ^ 

Y        =fYdt  v=/vdt 
o             -'      o  ^ 

r        =/rdt  \p      =      j   r  dt 

p        =/pdt  c|)      =      /pdt 

q        =/qdt  e=/qdt 
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A 


X 


NAV 


X 


NOTE 


Y 


1)  ip=0  when  vehicle  heads  parallel  to  +  X 
axis  . 


■NAV 


NAV 


2)  Vehicle  in  figure  above  shown  with  negative 
thrust  vector  angle  5,  positive  i/; ,  negative 
8  ,  negative  sway  velocity  v,  i.e. ,  in  a 
right  turn. 


Figure  1 
Definition  Of  Coordinate  System  (Part  I) 
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~/~  /-  g:  :1  ■' 


NOTE:   Direction  for  positive  velocity  and 
angular  rates  are  shown 

Figure  2 
Definition  Of  Coordinate  System  (Park  II) 
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BOW  ! 


^C^yVlv 


PORT 
SIDEWALL 


^Sx- 


^Sx-' 


w 


'  ^Cnvv!v[  P' 


'DY 


o 


Ha-r^ 


■*♦■ 


♦^s 


STARBOARD 
SIDEWALL 


r^io^i^^io 


TlO=°=*10 


N0TE:1)   p  and  p'  are  equivalent  point  force  centroids 

2)  5^,  5q  ,  6q,  5,  p,  are  thrust  vector  angles  and 

are  not  required  to  be  equal. 

3)  T„,  To,  T  ,  T  „  are  thrust  maenitudes  and  are 
not  required  To  be  equal. 

M-)   Force  directions  shown  are  for  a  right  turn 


Figure  3 
Surface  Effect  Ship  (Top  View) 


zzi^ 


T 
SIDE 

WATERLINE- 


xz   I 


NOTE 


1)   Right  Turn  Forces  Acting ,  i.e.,  5„,  S^,  5„, 
6  -,  p.  are  negative 

2'  "^SWE    '-   T7sin6y-^Tj,sin5g-^TgSin53^T^gSin6^g 


Figure  4A 
Surface  Effect  Ship  (.'-tern  View) 
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+  z 


•q 


FORW 


NOTE:       1)      TpQj^^=T^cos57-^TgCos6Q+TgCos6g+T^QCOs5^0 
(see    Fig.     3) 


Figure  UB 
Surface  Effect  Ship  (Side  View) 
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B.   FORCES  AND  MOMENTS 

1.   Surge  Forces  (See  Fig.  4B) 

The  forward  acceleration  was  determined  to  be 

generated  by  the  summation  of  thrust  vectored  in.  the  forward 

direction  T-^    .   The  resultant  acceleration  of  the  SES  was 
rorw 

assumed  to  be  opposed  by  a  retarding  force  exhibiting  a 
"velocity-squared"  characteristic.   Additionally,  Newton's 
laws  of  motion  specified  a  "v«r"  product  that  contributed 
to  the  drag  component  of  the  equation  when  the  vehicle  was 
in  a  turn  (centrifugal  force  reaction  term) .   The  simplified 
equation  of  motion  for  the  surge  acceleration  follows  from 
a  summation  of  forces  in  the  u  direction  as  depicted  in 
Fig.  4B. 


^  =   ^"^forw/"^^   "   ^^DX^^/^^   +  ^^ 


2.   Sway  Forces  (See  Fig.  4A) 

The  side  acceleration  of  the  vehicle  was  analyzed 
to  be  the  result  of  the  summation  of  the  thrust  vectored 
parallel  to  the  vehicle's  y-axis  and  a  retardation  "velocity 
squared"  phenomena  such  as  that  produced  by  the  cross  flow 
drag  term  of  the  sidewall  as  described  in  Ref.  6.   This 
retardation  force  is  augmented  by  a  "u-r"  product  as 
specified  by  Newton's  Law. 


^   =   ^side/"^   -   Sy  ^1^1/^   -   ^^ 
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2 

Note  the  requirement  to  force  the  C  v  /m  term  to 

maintain  the  sign  of  the  sway  velocity  such  that  both  left 

and  right  turns  may  be  computed. 

3.   Yaw  Moments  (See  Fig.  3) 

The  yaw  acceleration  was  a  turning  moment  summation 

which  was  defined  to  be  influenced  by  three  primary  forces 

operating  over  moment  arms.   It  is  of  interest  to  note  that 

the  third  term  is  an  "added  mass"  term  which  is  considered 

to  operate  on  the  same  pressure  point  as  the  C   v|v|  term. 

The  term  T    is  a  summation  of  moments  generated  by  the 
yaw  ^  -^ 

four  thrusters  as  shown  in  Fig.  3. 


-  ^2^10  ^°^  *10  "  ^o^7  ^^^  ^7  "  ^8  ^^^  ^8 

-  V9  ^^^  ^9  "  VlO  ^^^  '^lO 


^   =   Sy  ^  ^l^l/^Z   ^   Vw/'z   ^   ^22^^V^Z 


4.   Pitch  Moments  (See  Fig,  4B) 

The  pitch  acceleration  was  assumed  to  be  the  summa- 
tion of  moments  generated  by  forward  thrust,  ^^  ^^,  the 
buoyancy  of  the  sidewalls  of  the  SES ,  F^^  and  F^,    and  a 
vertical  force  generated  at  the  bow  of  the  vehicle,  F  . 
This  vertical  bow  force  was  defined  to  be  a  lumped  parameter 
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term  which  modeled  the  reaction  force  due  to  plenum  pressure 
acting  against  the  bow  seal. 


£ ,   =   average  draft  of  bow  seal 


^w,  =  II        +  ( (£   tan  e)/2)£ 
1     XX,       X       '    : 


A   =   £  £    -  ((£   tan  e)/2)£ 

W^       XX,         X         ^     X 


F,   =   A    £    pg 
1       w,   X2 


^2  =  ^2  "''^  P5 


^3  =  <=dpPb"e<^d  -  Si  *^"  «' 


£   =  average  draft  of  sidewall 
^1 


p,   =  plenum  pressure 


w,  =  average  wetted  sidewall  area  of  the  stern 
^x^  =  width  of  one  sidewall 


w   =  width  of  bow  seal 
e 


A   =  average  wetted  sidewall  area  of  the  bow 
^2 

A^^  =  added  mass  coefficient 


p   =  water  density 
g   =  32.3  ft/s^ 


£_^  =  lever  arm  of  bow  seal 
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S,~      =   lever  arm  of  buoyancy  force 

-  F^ij  -  A33Uq)/Iy 

Note  the  added  mass  term  A-^uq  which  was  required 
to  provide  damping  to  the  pitch  moment  and  is  specified 
in  Ref.  1. 

It  was  found  that  the  response  of  the  SES  to  a 
pitch  perturbation  without  the  added  mass  component  was 
undamped  and  approximately  sinusoidal.   Unlike  the  yaw 
equation  where  the  added  mass  term  was  a  small  contributor 
to  damping,  the  pitch  added  mass  term  was  found  to  be  of 
significant  importance  in  modeling  the  known  pitch  motion. 
5.   Roll  Moments  (See  Fig.  4A) 

The  roll  acceleration  equation  was  analyzed  to  be  a 

summation  of  moments  generated  by  buoyancy  forces  of  the 

port  and  starboard  sidewalls,  F   and  F  ,    the  thrust 
^  sp      ss 

vectored  parallel  to  the  y-axis  of  the  vehicle,  '^s±d&'    ^  side 

force  C_   v|vl,  and  a  lumped  coefficient  vertical  force, 

F   .   The  vertical  force  F   was  defined  as  the  force 
sw  sw 

generated  in  a  turn  due  to  the  sidewall  curvature  (dead 
rise  angle)  acting  against  the  cross  flow  of  water. 
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F     =   pg  A   £  , 
sp      ^   wp  dp 


F     =   pg  A    £ , 
ss      ^   ws   ds 


F    =   C    V  V 
sw      DZP   ' 


^j  =   average  draft  of  SES  sidewall 

il,  =      Ij      -    i        tan  d) 
dp      d,     sw      ^ 

£,  =   £,   +  £    tan  d) 
ds      d,     sw     ^ 


A    =   average  wetted  area  port 


A    =   average  wetted  area  starboard 
ws  ^ 


p=   ((F    -   F      )l        -T.,   £+  C^,,  vlvl  i 
^  sp     ss   sw     side   p     DY   '  '   2 


-  F    i         -   A-T  up)/I^ 
sw   sw     31   '^  '  X 

Again  note  the  added  mass  term  A--up/I  which  was 
found  to  be  essential  in  modeling  the  damping  phenomena  of 
the  roll  motion. 

In  summary,  the  simplified  5  DOF  equations  of  motion 
used  in  RTS5D  for  the  3K  TON  SES  are: 


•  2 

SURGE     u   =   T-    /m  -  C^^u  /m  +  vr     I Fig.  4bJ 
forw      DX  ' 


SWAY      V  =   Tg^^^/m  -  Cj^^  v  |  v  | /m  -  ur     iFig.  4AJ 
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YAW 


^      =      Vw/'Z    -^    ^DY    ^1^1     V^Z    -^   ^22^^V^Z  f^^^-    3^ 


-   F^Z^    -   A^^uq)/!^  [Fig. .4BJ 


ROLL  p=       ((F         -F)£         -T.,       il+    C^„   vIvU 
^  '  '    sp  ss      sw  side     p  DY      '     '    z 


-    ^sw^sw   -   ^3    ^P^/^X  f^^^-    ^^1 


D.   PARAMETER  IDENTIFICATION 

These  equations  are  an  extension  of  the  3  DOF  flat  turn 

SES  model  developed  by  Gerba  and  Thaler  in  Ref.  4.   The 

identification  of  craft  parameters  C_„,  C_-, ,  and    is 

Da    DY       w 

described  in  Ref.  4  and  repeated  here  for  completeness. 

The  surge  drag  coefficient  C   is  determined  by  selecting 

DX 

a  steady  state  turn  condition,  where  T^    ,  u,  v,  r,  and  m 
are  known. 

u  =   0   =   -^2£^  -  -^^^—     +  vr 
m         m 


from  which 


T^    +  mvr 

^DX  "  2 

u 


29 


The  sway  drag  coefficient  C   is  determined  by  the  s 


DY 


ame 


steady  state  turn  condition  which  requires 


where 


T  .  ^      C 
n               side      DY   I  I 
V  =   0   =  -   V  V   -  ur 


<=DY 


T  .  ,   -  mur 
side 


v|v| 


The  sway  drag  moment  arm  follows  utilizing  the  yaw 
acceleration  equation 

'^.r^r.        C_^  v|v|£    A_-  uv  Jl 
r  =   0   =   ^        +  '  '  w  ^22     w 


h  H  h 


which  yields 


-T 
Zaw. 


^       DY  v|v|  +  A22^v 


The  addition  of  the  pitch  and  roll  equations  introduces 

additional  lumped  coefficients  C   and  C   _.   These  are 

easily  solved  using  known  constants  1,1,    C^„,  Z    ,    i      , 
^  ^  w    p    DY    z    sw 

^^/  Pw/  ^j  ^nd  steady  state  values  of  F   ,  F  „,  F,,  F_, 
e   D   d  -^  sp   ss    1    2 

T-    ,  T  . ,  ,  u,  V,  and  r.   It  is  significant  to  note  that 
forw   side  ^ 

in  a  steady  state  condition  the  angular  roll  rate,  p,  and 
angular  pitch  rate,  q,  are  both  required  to  equal  zero. 
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This  is  in  contrast  to  the  yaw  rate,  r,   where  a  finite 
steady  state  value  is  desired  in  a  turn.   Thus  the  added 
mass  terms  of  A^^uq  and  A_-up  in  the  pitch  and  roll  equations 
are  not  utilized  in  the  determination  of  craft  parameters; 
their  function  is  strictly  confined  to  damping  of  their 
respective  accelerations.   Therefore,  utilizing  the  pitch 
equation  with  known  steady  state  turn  values. 


q   =   0   =   (T^     I      +  C^„p,w  (jl,  -  £-,  tan  9)  £_. 
^  forw   p    DP^b  e  d    31         31 


^  ^2^  -  Sx^^^2  -  V3^/^Y 


yields 


Sp 


2 

Sx^  ^z  "^^1^3  "  -^2^3  '  "^forw  ^p 

p,w  {l^   -  £_,  tan  e)Ji-,, 
^b  e   d    31        31 


The  roll  equation  under  steady  state  turn  conditions 
yields 


from  which 


SZP  =   ((^sp  -  ^ss>^w  -  ^side  V  ^  Sy  ^ssl^ssl^z'/ 


(v   Iv  \     I   ,) 
^    ss '  ss '   sw' 
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where 


F  =  port  sidewall  buoyancy  in  pounds     [Fig.  4A] 
sp 

F  =  starboard  sidewall  buoyancy  in  pounds   [Fig.  4A] 

T  .  ,  =  thrust  vectored  parallel  to  vehicles  y-axis 
in  pounds   [Fig.  3] 

V  =  sway  steady  state  velocity  in  feet/sec 

C_^y  =  sway  drag  coefficient  for  steady  state  condition 

Z  =  sway  drag  moment  arm  in  feet    [Fig.  3] 

i  =   effector  thrust  moment  arm  in  feet   [Fig.  4A] 
P 

£  =  surge  drag  moment  arm  in  feet    [Fig.  4A] 
z 


sw 


=  sidewall  buoyancy  force  moment  arm  in  feet 
[Fig.  4A] 


32 


III.   INTEGRATION  METHOD 

A.   ERROR  ANALYSIS 

A  rectangular  (RECT)  type  integration  was  used  to 
compute  the  velocity  and  positional  variables  of  the  force 
equations.   Selection  of  the  type  of  integration  over  a 
more  complex  type  such  as  Runga  Kutta,  Variable  Step  enabled 
the  loop  delay  time  due  to  integration  to  be  held  to  a 
minimum . 

An  analysis  of  the  error  introduced  into  the  solution 
of  the  simplified  5  DOF  equations  through  the  use  of 
rectangular  integration  vice  Runga-Kutta,  Variable  Step  was 
made.   The  simplified  5  DOF  equations  of  motion  were 
programmed  on  the  NPS  IBM  36  0/67  computer  using  the 
Continuous  System  Modeling  Program  (CSMP) .   One  series  of 
computer  runs  was  made  utilizing  rectangular  integration 
over  a  spectrum  of  step  sizes  ranging  from  0.01  sees  to 
0.25  sees.   A  second  computer  run  was  made  utilizing  Runga- 
Kutta,  Variable  Step  type  integration.   The  steady  state 
results  are  shown  in  Table  I.   Examination  of  the  final 
value  magnitudes  established  that  no  significant  steady 
state  error  was  due  to  rectangular  integration  for  the  time 
increments  utilized. 

Transient  error  of  the  variables,  r,  (t> ,    and  6  was 
examined.   A  series  of  computer  runs  using  incrementally 
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larger  rectangular  integration  step  sizes  was  compared 
against  an  identical  run  using  Runga-Kutta  Variable  Step 
type  integration.   The  first  peak  absolute  magnitude 
results  are  shown  in  Table  II.   The  RTS5D  had  an  iteration 
time  of  .084  seconds  therefore  the  "transient"  magnitudes 
were  in  error  by  approximately  2%  to  3%  due  to  rectangular 
integration. 
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TABLE  II.   TRANSIENT  FIRST  PEAK  ERROR  ANALYSIS 
FOR  RECTANGULAR  INTEGRATION 


DELT           r                   (J)  0 

1.64%  2.58% 

2.41%  2.84% 

3.38%  3.25% 

3.65%  8.57% 


.050 

1.42% 

.075 

1.89% 

.100 

2.16% 

.125 

2.35% 

Note:     1)   Step  effector  angle  input  of  25° 
and  velocity  (u)  of  56  knts. 

2)   All  percentages  are  referenced  to 

Runga-Kutta  variable  step  magnitudes. 
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IV.   IMPLEMENTATION 

A.   PHILOSOPHY 

Given  the  decision  to  develop  a  real  time  motion  analysis 
simulator  for  the  3K-SES,  the  various  options  available 
were  examined.   Prime  requirements  for  such  a  simulator 
were  the  necessity  for  real  clock  interfacing,  suitable 
input-output  ports,  and  a  dedicated  Central  Processing 
Unit  (CPU) . 

The  "all  analog"  computer  simulation  was  particularly 
attractive  from  the  aspect  of  interfacing  man  generated 
control  inputs  into  an  on-going  real  time  solution.   With 
the  analog  approach  instantaneous  solutions  are  obtained 
from  patchboard  amplifiers.   Two  considerations  which  did  not 
lend  themselves  well  to  the  all  analog  technique  were  how  to 
maintain  high  accuracy  solutions  in  wide  ranging  variables 
when  constrained  to  a  voltage  limit  of  plus  or  minus  100 
volts,  and  how  to  output  the  system  variables  in  such  a 
manner  that  a  man  could  physically  interpret  what  was 
happening  and  respond  with  the  required  corrective,  timely 
control  effort.   Additional  supplementary  problems  associ- 
ated with  the  "quarter-square"  multiplication  error  and 
amplifier  drift/run  duration  limitations  gave  rise  to 
examining  other  possible  approaches. 

The  "all  digital"  computer  simulation  offered  high 
accuracy  with  wide  ranging  variables  and  a  capability  to 
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interface  with  two  graphics  display  consoles.   Additionally 
a  line  printer  was  available  to  provide  "hard  copy"  of 
variables  of  interest  as  the  solution  was  being  processed. 
Two  major  disadvantages  were  the  serial  and  discrete  nature 
of  the  integration  method  of  a  digital  computer,  and  the 
lack  of  high  voltage,  high  accuracy  potentiometer  controls 
to  allow  the  man  to  influence  the  motion  of  the  vehicle  as 
the  solution  was  being  presented  to  him. 

A  desire  to  take  advantage  of  the  best  features  of  both 
the  analog  and  digital  alternatives  led  to  a  decision  to 
develop  the  real  time  simulator  (RTS5D)  of  the  3K-SES 
utilizing  the  simplified  5  DOF  equations  of  motion  developed 
in  Section  II. 

B.   REQUIREMENTS 

The  following  criteria  were  required  of  the  RTS5D: 

1)  The  simplified  equations  of  motion  would  be  solved 
and  the  results  output  on  a  real  time  basis. 

2)  The  solution  would  be  subject  to  real  time  control 
efforts  generated  by  a  man  observing  the  output. 

Development  of  the  simplified  5  DOF  equations  of  motion 
and  the  use  of  rectangular  type  integration  for  simulating 
the  3K-SES  ensured  that  loop  delay  attributable  to  digital 
computational  needs  would  be  minimal  and  thus  the  solution 
would  be  compatible  to  a  real  time  environment. 

The  requirement  to  output  the  results  on  a  real  time 
basis  such  that  a  man  observed  and  responded  with  real  time 
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control  efforts  necessitated  the  use  of  graphics  terminals 
interfaced  with  a  digital  computer.   Large  volumes  of  data 
and  graphical  cues  were  presented  to  the  pilot  and  updated 
on  a  real  time  basis  in  much  the  same  manner  as  found  on 
the  bridge  consoles  of  the  3K-SES  shown  in  Ref,  7.   Software 
algorithms  designed  to  present  "alpha-numeric"  data  and 
graphic  perspectives  corresponding  to  the  vehicles  current 
state  variables  enabled  the  observer  to  respond  on  a  real 
time  basis. 

Real  time  man  generated  control  requirements  necessitated 
the  need  for  a  hybrid  environment  and  its  associated  inter- 
face between  a  digital  computer  and  an  analog  computer. 
The  RTS5D  utilizes  four  (4)  computer  systems:   One  (1) 
medium  size  general  purpose  digital  computer  (XDS-9300) , 
two  (.2)  small  graphics  digital  computers  (AGT-10)  ,  and  one 
(1)  analog  computer  (CI5000) .   The  XDS-9300  has  a  CPU 
capable  of  handling  32K  words  of  main  (core)  memory  with  a 
magnetic  drum  for  secondary  storage.   The  XDS-930  0  may  be 
programmed  in  assembly  language  or  FORTRAN.   The  CI-5000 
is  an  analog  computer  with  a  full  complement  of  digital 
logic.   Hardwired  trunk  lines  linking  the  XDS-9300  and 
CI-5000  allow  hybrid  computer  operations.   The  two  AGT-10 's 
are  small  general  purpose  digital  computers  with  8K  of  main 
memory  and  magnetic  disk  for  secondary  storage.   The  AGT-10 
systems  have  the  capability  of  accepting  input  and  of 
displaying  output  on  their  CRT  screens.   Both  the  CI-50  00 
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and  the  AGT-10  have  the  capability  of  communicating  with 
the  XDS-93  00  directly,  but  cannot  communicate  with  each 
other  except  via  the  XDS-9300  [Ref.  8]. 

C.  HARDWARE  DESCRIPTION 

Hardware  layout  for  the  RTS5D  is  shown  in  Figure  5A. 
A  "joystick"  type  high  voltage  potentiometer  attached  to 
the  pilot's  seat  was  hardwired  to  the  CI-50  00.   This  joy- 
stick served  the  purpose  of  the  pilot's  effector  angle 
controller.   Left  or  right  stick  transmitted  an  appropriate 
nozzle  deflection  angle  in  the  form  of  a  voltage  level  to 
the  CI-5000  where  it  was  quantized  in  an  analog  to  digital 
process  and  introduced  into  the  XDS-9300 's  ongoing  solution 
of  the  simplified  5  DOF  equations.   In  like  manner  the 
voltage  levels  of  six  other  "dial"  potentiometers  used  to 
control  thruster  output  levels  and  program  modes  were 
introduced  into  the  solution  process  as  deemed  appropriate 
by  the  operator.   The  six  potentiometers  were  located  on  a 
"thruster  console"  (See  Fig.  5B)  within  reach  of  the  pilot 
{See  Fig.  5A) .   Additional  hardware  included  a  standard 
digital  computer  line  printer,  an  analog  8  channel  strip 
chart  recorder,  and  an  analog  X-Y  plotter, 

D.  SOFTWARE  DESCRIPTION 
1.   Introduction 

The  RTS5D  software  package  was  a  FORTRAN  IV  digital 
computer  program  which  consisted  of  a  main  program  and  two 
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Figure  5A 
RTS5D  Block  Diagram 
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Figure  SB 
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major  subroutines.   The  program  gave  the  operator  the  option 
to  start,  stop,  and  hold  the  dynamics  solution  process. 
The  software  package  controlled  a  multiplexing  algorithm 
used  on  the  graphics  display  consoles  and  sampled  all 
analog  control  and  timing  inputs.   A  program  flow  chart  is 
shown  in  Fig.  6A  and  a  detailed  flow  chart  for  the  graphics 
multiplexing  is  shown  in  Fig.  6B. 

The  program  began  by  initializing  all  program 
variables  and  parameters  and  then  requested  from  the  operator 
any  parameter  changes.   The  request  and  input  of  changes 
were  made  on  a  TTY-35  Teletypewriter  (See  Fig.  5A) .   The 
program  then  started  the  real  time  1000  Hz  clock  timing 
circuits  in  the  CI-5  000  analog  computer  and  entered  the  main 
program  iteration  loop.   In  the  main  loop  the  sampling  of 
analog  control  inputs  by  the  pilot  was  accomplished  by 
subroutine  ADDA.   Craft  force  equations  were  computed  and 
integrated.   A  perspective  view  of  the  simplified  form  of 
the  SES  and  an  electronic  road  were  utilized  in  the  RTS5D 
to  give   real  time  visual  cues  to  the  pilot  on  the  status 
of  the  vehicle's  states.   This  section  of  the  display  gave 
the  pilot  his  primary  guidance  information  and  was  updated 
on  the  pilot's  graphics  display  console  on  every  iteration 
of  the  main  program  loop.   The  program  then  evaluated 
whether  an  accidental  failure  mode  had  been  selected  by  the 
operator.   If  the  failure  mode  was  selected  the  program 
searched  the  craft's  states  to  see  if  the  recovery  criteria 
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Figure  6A 
RTS5D  Program  Flow  Chart 
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Figure  6B 
RTS5D  Graphics  Display  Multiplex  Algorithm 
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had  been  satisfied  through  the  maneuvering  inputs  of  the 
pilot.   If  the  recovery  criteria  was  fulfilled,  the  program 
went  into  a  hold  status  and  awaited  further  instructions 
by  the  operator  on  when  to  restart  the  run.   If  the  recovery 
criteria  had  not  been  fulfilled  the  program  continued  the 
main  loop  by  incrementing  a  graphics  multiplex  index  and 
entering  the  graphics  display  multiplex  algorithm  to  search 
for  and  display  the  appropriate  information  array  associated 
with  the  current  value  of  the  graphics  index  I  (See  Fig.  6B)  , 
A  graphics  multiplexing  algorithm  was  established  to  reduce 
the  iteration  time  required  for  the  main  loop.   Each  informa- 
tion array  contained  approximately  50  elements  associated 
with  positions  on  the  display  screens.   Each  array  had  a 
display  time  which  was  required  to  encode  and  project  the 
array  on  the  screen.   This  display  time  varied  directly 
with  the  number  of  words  in  the  array.   The  prioritization 
of  the  graphic  display  console  information  arrays  was  made 
such  that  the  perspective  of  the  SES  and  the  electronic 
road  were  displayed  on  each  iteration  and  were  not  subject 
to  the  program's  graphics  display  multiplexing  routine. 
All  other  information  arrays  were  multiplexed  on  a  time 
share  basis  such  that  on  any  main  loop  only  two  information 
arrays  were  projected  on  the  screen,  the  SES/road  and  the 
array  designated  by  the  graphics  display  index  I.   This 
concept  was  mandatory  due  to  the  large  volume  of  data/ 
graphics  being  presented.   It  is  instructional  to  recall 
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the  fact  that  all  variables  were  computed  and  updated  in 
the  memory  of  the  XDS-9300  on  every  loop,  whereas  the 
display  routine  was  multiplexed  on  a  time  share  basis. 
The  main  loop  iteration  time  was  84  msec.   Information 
arrays  in  the  multiplex  routine  were  updated  on  the  graphics 
display  consoles  every  2.0  sec.   On  exit  from  the  graphics 
display  multiplex  algorithm  the  program  read  the  value  of 
the  real  time  clock  and  then  reset  the  clock  to  zero.   From 
this  the  A  time  magnitude  was  determined.   Following  this 
the  program  searches  for  print,  restart,  and  hold  flags 
which  are  set  by  the  operator.   The  program  loops  back  to 
the  AD/DA  subroutine  and  begins  a  new  iteration. 
2 .   Graphics  Techniques 

The  two  AGT-10  Graphics  Display  Consoles  were 
utilized  to  display  real  time  data  to  the  operator  such  that 
he  could  control  the  craft.   The  process  of  transmitting  and 
displaying  desired  data  to  an  AGT-10  was  accomplished 
utilizing  FORTRAN  callable  graphics  subroutines  in  the 
XDS-9300.   One  graphics  subroutine  TEXTO  was  used  to  transmit 
alpha-numeric  data  and  one  graphics  subroutine  GRAPHO  was 
used  to  transmit  graphical  data. 

Subroutine  TEXTO  transmitted  alpha-numeric  data  to 
the  AGT-10  screen.   The  positioning  of  the  desired  data  on 
the  screen  itself  was  accomplished  by  utilizing  a  TEXTO 
screen  coordinate  system  of  dividing  the  12  inch  by  12  inch 
AGT-10  screen  into  40  equally  spaced  lines  from  top  to 
bottom  with  each  line  having  room  for  96  individual 
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characters  on  it  as  shown  in  Fig.  7A.   Subroutine  TEXTO 
statements  analogous  to  FORTRAN  "WRITE"  and  "FORMAT"  state- 
ments were  used  to  specify  the  letter  character  or  variable 
magnitude  to  be  displayed  and  their  position  on  the  display 
screen. 

Subroutine  GRAPHO  was  used  to  transmit  the  "picture" 
portion  of  the  output  to  the  display  screen.   The  coordinate 
system  utilized  by  GRAPHO  is  a  normalized  X-Y  system  with 
the  ordinate  and  abscissa  intersecting  in  the  center  of  the 
screen  as  shown  in  Fig.  7B.   A  point  on  the  GRAPHO  coordinate 
system  corresponding  to  (l.,0.)  would  be  located  at  the  right 
edge  of  the  screen  centered  between  the  top  and  bottom  of 
the  screen.   A  line  could  be  generated  on  the  screen  by 
specifying  the  line  start  coordinates,  line  stop  coordinates 
and  commanding  a  "draw".   The  AGT-10  would  connect  the  two 
designated  points  and  a  line  would  be  displayed  on  the  screen. 
Concurrently  if  only  two  point  sources  of  light  had  been 
desired  to  be  displayed,  the  command  "jump"  could  be 
utilized  such  that  the  AGT-10  would  not  connect  the  second 
coordinate  point  with  the  first  coordinate  point.   In  this 
manner  a  series  of  commands  was  generated  via  subroutine 
GRAPHO  whereby  a  picture  was  drawn  on  the  AGT-10  screen 
using  an  array  of  sequential  coordinate  points  and  their 
respective  "draw"  or  "jump"  instruction. 
3.   Analyst  Graphic  Display  Console 

The  display  presented  to  the  analyst  during  the 
conduct  of  a  run  is  shown  in  Fig.  8,   The  magnitudes  of  the 
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Figure  7A 
Subroutine  Texto  Coordinate  System 
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Figure  8 
Analyst  Graphic  Display  Console 
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the  variables  were  computed  by  the  XDS-9300  on  every  iteration 
of  the  program.   The  XDS-9300  then  followed  a  display  multi- 
plexing algorithm  (See  Fig.  6B)  which  updated  the  variables 
every  2.  seconds.   The  variables  were  categorized  by  position, 
velocity,  and  acceleration.   Each  category  had  three  lines 
of  data  displayed.   The  "present"  data  line  displayed  the 
real  time  magnitude  of  the  variable  column  it  was  located 
in.   The  "maximum"  and  "minimum"  were  the  max  and  min 
magnitudes  that  had  been  computed  for  the  variable  during 
that  particular  run.   At  the  start  of  a  new  run  all  variables 
are  initialized  (See  Fig.  6A) . 

4.   Pilot  Graphics  Display  Console 

The  display  presented  to  the  pilot  during  the  conduct 
of  a  run  is  shown  in  Fig.  9A  and  Fig.  9B.   The  pilot  display 
contained  three  major  sections,  VARIABLE  DYNAMICS,  SEA  TRACK/ 
EMERGENCY  PERFORMANCE  INDEX,  CRAFT  PERSPECTIVE/NORTH  ROAD/ 
NAV  DATA. 

a.   Variable  Dynamics 

The  time  histories  of  the  accelerations  and  their 
respective  velocities  were  presented  in  graphical  form.   An 
acceleration  plot  and  its  respective  velocity  plot  shared 
the  same  electronic  strip  chart  such  that  one  was  super- 
imposed over  the  other.   The  acceleration  plots  were 
displayed  as  "light"  toned  lines  whereas  the  velocity  plots 
were  displayed  as  "heavy"  toned  lines.   As  the  RTS5D  program 
incremented  in  real  time,  it  presented  an  incremental 
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Figure  9-B 
Pilot  Graphic  Display 
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increase  in  each  displayed  state  time  history.  The  time 
histories  would  progress  in  a  left-to-right  movement  and 
upon  reaching  the  right  most  limit  of  the  VARIABLE  DYNAMICS 
section  were  erased  from  the  screen  and  restarted  at  the 
left  side  of  the  electronic  strip  chart.  Execution  of  a 
maneuver  of  the  vehicle  would  immediately  be  displayed  in 
the  time  histories  displays. 

b.   Sea  Track/Emergency  Performance  Index 

The  sea  track  portion  of  this  display  area  was 
a  navigational  plot  of  the  ship  trajectory  (X  ,Y  )  (See 
Fig.  9B) .  This  track  was  updated  every  2  seconds  as  imple- 
mented in  the  XDS-9300  multiplex  algorithm.   The  track 
would  erase  itself  and  begin  anew  at  the  craft's  current 
position  every  3.5  minutes. 

The  performance  of  the  pilot  during  an  acci- 
dental failure  situation  was  measured  and  graded  through 
the  use  of  three  performance  indices  J.  (TRACK  DEVIATION) , 
J2  (CONTROL  EFFORT) ,  and  J^  (OVERALL  PERFORMANCE) .   All 
performance  tests  were  initialized  such  that  the  vehicle 
was  on  a  heading  of  000°  (NORTH)  and  on  the  X^-^„  axis  with 

NAY 

a  steady  state  velocity. 

J, ,  TRACK  DEVIATION,  was  an  "integral-error- 
squared-plus-time"  performance  index  which  penalized  the 
distance  that  the  pilot  allowed  the  vehicle's  e.g.  to 
drift  away  from  the  X    axis  in  addition  to  the  length 
of  time  the  vehicle  remained  off  the  X    axis  (See  Fig.  10) 
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Figure  10 
Track  Deviation  Perfarmance  IndexCJ  ) 
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At  the  systems  entry  into  an  accidental  failure  mode  of 
operation  the  XDS-9300  computed  J^  as 

final      2 

Jt   =    /        (Y  )   dt   +   TIME 
1      ^^  o 


J 21    CONTROL  EFFORT,  was  an  index  used  to 
penalize  the  pilot  for  excessive  control  effort  during  the 
man  controlled  recovery  maneuver  segment  of  the  accidental 
failure  test.   The  J2  index  was  an  integral  "error-squared" 
index,  as 

t, 


'final 

I 

t 


J,   =    /        (6)^  dt 


Note  that  the  6  under  the  integral  sign  is  the  6  associated 
with  the  nozzle  deflection  of  the  three  operational  thrusters 
of  the  SES,  and  does  not  include  the  nozzle  angle  of  the 
failed  thruster  (i.e.,  starboard  thruster) .   A  graphical 
interpretation  of  J^  is  shown  in  Fig.  11. 

J^,  OVERALL  PERFORMANCE,  was  a  magnitude 
representing  the  sum  of  J,  and  J^ 


J3   =   J^  +  J2 
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Figure  11 
Control  Effort  Performance  Index(J„) 
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At  the  completion  of  the  test  the  performance  algorithm 
in  the  XDS-9300  was  programmed  to  display  the  numeric 
scores  representing  J,  and  J-  and  the  word  evaluation, 
EXCELLENT,  ABOVE- AVERAGE ,  or  GOOD  below  the  overall  performance 
heading  shown  in  Fig.  9B.   The  magnitude  range  of  the 
performance  index  J^  for  each  of  the  3  word  evaluations  are 
shown  in  Table  III. 


TABLE  III.   OVERALL  PERFORMANCE 
INDEX  RATINGS 


RATING        EXCELLENT     ABOVE  AVERAGE        GOOD 
RANGE         J^  <  250      250  1  J3  <  500       500  <_  J 


"c.   Craft  Perspective/North  Road/NAV  Data 

Real  time  visual  cues  to  the  pilot  were  generated 
on  the  lower  half  of  the  pilot's  graphics  display  console 
in  the  form  of  a  perspective  view  of  a  simplified  shape  of 
the  craft  from  the  rear  looking  down  and  forward  onto  the 
craft  and  the  area  immediately  ahead  of  the  craft's  path. 
An  electronic  north-heading  road  shown  forward  of  the  craft 
provided  information  to  assist  the  operator  in  proper 
control  of  the  vehicle.   Numeric  navigation  data  concerning 
the  craft's  position  and  the  status  of  the  nozzle  angles 
and  effector  thrust  settings  were  presented  in  the  area 
immediately  below  the  craft  perspective  (See  Fig.  9A  and 
Fig.  9B) . 
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The  craft  perspective  view  presented  real  time 

graphical  information  on  the  vehicle's  turning  rate  (r) , 

velocity  (u)  ,  pitch  angle  (9),  and  roll  angle  (({)).   The 

electronic  north  road's  centerline  was  the  positive  X,,,,, 

^         NAV 

axis   (See  Fig.  1) .   The  north  road  was  displayed  on  the 

pilot's  screen  whenever  the  pilot  had  maneuvered  the  craft 

to  within  1500  feet  of  the  X„,..  axis  and  the  heading  of  the 

NAV  ^ 

vehicle  was  on  a  northerly  course.   At  other  times  the  north 
road  was  not  visible  to  the  pilot.   When  displayed,  the 
north  road  provided  the  pilot  with  visual  cues  on  the 
vehicle's  heading  and  navigation  position  with  respect  to 

the  fixed  X„,„  axis.   When  the  pilot  maneuvered  the  craft 

NAV  ^ 

such  that  the  north  road  was  not  visible  graphical  turn  rate 
information  was  determined  from  a  zero  turn  rate  reference 
line  located  just  above  the  bow  of  the  SES.   The  pilot 
established  a  zero  turn  rate  whenever  the  center  point  of 
the  craft's  bow  was  lined  up  just  below  the  zero  turn  rate 
reference  line.   Variable  speed  lines  moving  from  the 
craft's  bow  toward  the  stern  gave  a  graphic  representation 
of  the  forward  velocity  u  (See  Fig.  12) .   A  sequence  of 
displays  during  a  pilot  commanded  control  effort  of  left 
nozzle  deflection  and  then  right  nozzle  deflection  is  shown 
in  Fig.  13,  Fig.  14,  and  Fig.  15. 

The  north  road  positioning  on  the  pilot's 
graphic  display  console  was  solved  using  the  following 
assumptions : 
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NOTE:   1)  Craft  "shows  center  of  bow  is  left  of 

reference  line  hence  a  left  turn  status. 


Figure  12 

SES  Perspective  Information 
Without  North  Road 
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VARIABLE    DYNAMICS 


U         (HEAVY) 
UDOTCLIGHT) 

V         (HEAVY) 
VDOT( LIGHT) 

R         (HEAVY) 
RDOT( LIGHT) 

P         (HEAVY) 
PDOT(LIGHT) 

Q         (HEAVY) 
QDOT(LIGHT) 


PL 


^^k. 


SEA  TRACK 
I 


-^ 


EMERGENCY    PERFORMANCE    INDEX 
TRACK         CONTROL  OVERALL 

DEVIATION      EFFORT    PERFORMANCE 


CONTROL    INPUTS  TIME 

THRUST         RUDDER  DELT   ELAPSED 

400000  .32  .084  10 

THRUSTl    THRUST2    THRUSTS    THRUST4 
100000      100000      100000      100000 


NAVIGATIONAL    DATA 

PITCH   HEADING/RATE    DRIFT   ROLL   SPEED 

.52  355/-1.0         .02  -.02  56 

RUDl  RUD2  RUD3  RUD4 

.32  .32  .32  .32 


NOTE:  1)  Craft  has  been  commanded  by  pilot  into  a  left  turn 
by  deflecting  nozzles  .32  radians  (+18°).   Vehicle 
is  entering  left  turn  and  is  passing  through  a  heading 
of  355°  with  a  turn  rate  of  -l°/sec  (see  NAVIGATION 
DATA  readout) . 

2)  SEA  TRACK  and  electronic  road  show  vehicle  has  dis- 
placed slightly  left  of  north  road  centerline . 

3)  VARIABLE  DYNAMICS  show  velocity  and  acceleration  time 
histories  being  generated  as  manuever  develops . 


Figure  13 
SES  MANUEVER  (I) 
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VARIABLES 


U        (HEAVY) 
UDOT(LIGHT) 

V         (HEAVY) 
VDOT( LIGHT) 

R         (HEAVY) 
RDOT( LIGHT) 

P         (HEAVY) 
PDOT( LIGHT) 

Q         (HEAVY) 
QDOT(LIGHT) 


£^ 


SEA  TRACK 
I 


^ 


R^ 


itX*. 


EMERGENCY   PERFORMANCE    INDEX 
TRACK         CONTROL  OVERALL 

DEVIATION      EFFORT    PERFORMANCE 


CONTROL    INPUTS  TIME 

THRUST         RUDDER  DELT   ELAPSED 

400000  .32  .084  20.2 

THRUSTl    THRUST2    THRUST 3   THRUST4 
100000      100000      100000      100000 


NAVIGATIONAL    DATA 

PITCH   HEADING/RATE    DRIFT   ROLL  SPEED 

.58                320/1,2         .04      -.5  51 

RUDl              RUD2                 RUD3  RUD4 

.32                .32                  .32  .32 


NOTE:  1)  Pilot  has  held  left  nozzle  deflection  angle  of  .32 
rad  (+18°)  and  vehicle  has  displaced  left  of  north 
road.   Vehicle  has  heading  of  320°  or  40°  off  the 
north  road  centerline. 

2)  SEA  TRACK  confirms  trajectory  of  craft  in  left  turn 

3)  VARIABLE  DYNAMICS  show  +u^+v,-r,  and  q  in  steady 
state  while  roll  rate  p  is  still  seeking  steady 
state  condition. 

Figure  14 
SES  Manuever  (II) 
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U    (HEAVY) 
UDOT( LIGHT) 

V    (HEAVY) 
VDOT( LIGHT) 

R    (HEAVY) 
RDOT( LIGHT) 

P    (HEAVY) 
PDOT( LIGHT) 

Q    (HEAVY) 
QDOT(LIGHT) 


VARIABLE  DYNAMICS 

^1  I 


SEA  TRACK 


^17 


Co^^K\o^ 


EMERGENCY   PERFORMANCE    INDEX 
TRACK         CONTROL  OVERALL 

DEVIATION      EFFORT    PERFORMANCE 


CONTROL    INPUTS  TIME 

THRUST         RUDDER  DELT   ELAPSED 

400000  -.5  .084  33. 

THRUSTl    THRUST2    THRUST3    THRUST4 
100000 


NAVIGATION   DATA 

PITCH   HEADING/RATE    DRIFT    ROLL  SPEED 

.51               030/1.8        -.3         .8  50 

RUDl       RUD2      RUD3  RUD4 

-.5       -.5      -.5  -.5 


NOTE:  1)  Vehicle  has  been  manuevered  into  right  turn. 

2)  Vehicle  heading  0  30°  on  path  which  will  intersect 
north  road. 


Figure  15 
SES  Manuever  (III) 
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1)  The  view  of  the  pilot  was  restricted  to  plus  and 
minus  60  degrees  of  the  heading  of  the  craft.   This 
corresponded  to  the  far  right  and  far  left  edges 

of  the  pilot's  display  screen,  respectively. 

2)  The  starting  edge  of  the  road  would  remain  fixed  at 

a  point  ahead  of  the  vehicle  whenever  the  vehicle  was 
positioned  such  that  the  pilot  could  see  the  road  in 
his  restricted  viewing  arc. 

3)  The  road  would  give  visual  cues  on  the  craft's  left- 
right  displacement  from  the  road  centerline  and  on 
the  relationship  between  the  north  heading  of  the 
road  and  the  heading  of  the  vehicle. 

4)  The  road's  shape  as  it  appeared  to  the  pilot  would  be 

a  wide  path  terminating  at  an  apex  point  in  the  distance 
A  non  linear  relationship  was  used  to  compute  the 
position  of  the  north  road  (See  Fig.  16A  and  Fig.  16B) .   The 
Y  ^  displacement  of  points  A,  B,  and  C  on  the  graphics 
display  console  was  made  a  function  of  Y   only.   The  Y  ^ 
displacement  of  point  D  on  the  graphics  display  console  was 
made  a  function  of  Y   and  the  arc  tangent  of  the  craft's 
heading.   The  X  ^  displacement  of  the  points  A,  B,  C,  and 
D  were  fixed  at  a  constant  value,  thus  recalling  Fig.  7B, 
the  computation  of  the  coordinates  of  the  points  A,  B,  C, 
and  D  was 
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X 


NAV 


FAR  DISTANT  ROAD  APEX 


X 


o 


START  OF  ROAD 


VIEWING  ARC 
OF  PILOT 


^^NAV 


Figure  16A 
North  Road(I) 
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^NAV   ^^^^ 


•NAV 


Figure    16B 
North   Road (II) 
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^A   =   -(^o  ^  ^1^  ^A=  ^-Ol 


^B   =   -^o  ^B  =  ^-01 


Xj^   =   -(Yq  +  10000  (tan"-"-  i|;)  )  Y^^  =  +.1 


The  displacement  of  points  A,  B,  and  C  due  to 
vehicle  heading  were  neglected. 

E.   OPERATIONAL  MODES 

The  RTS5D  operated  in  any  combination  of  modes  as 
designated  by  the  operator  (See  Fig.  17) .   Mode  selection 
was  accomplished  via  "on-off"  switches  located  at  the 
console  of  the  XDS-9300  and  at  the  thruster  console  located 
at  the  pilot's  seat  (See  Fig.  5A  and  Fig.  5B) .   The  RTS5D 
ran  in  any  combination  of  the  modes  available  with  the 
"RESTART"  mode  taking  priority  over  all  other  modes  selected, 

The  FIXED  THRUST  and  FIXED  EFFECTOR  ANGLE  removed  the 
appropriate  man-generated  control  and  utilized  a  constant 
value  for  T  and  6  respectively. 

The  PRINT  mode  allowed  variables  to  be  printed  on  the 
XDS-9300  line  printer  at  intervals  selected  by  the  operator. 

The  HOLD  mode  put  the  program  into  a  loop  whereby  the 
display  screens  would  continue  to  present  the  last  image 
computed  until  such  time  as  the  mode  was  deselected.   At 
this  time  the  program  would  continue. 
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MODE  DESCRIPTION 

1  THRUST  MAGNITUDE  FIXED 

2  THRUST  MAGNITUDE  VARIABLE 

3  EFFECTOR  DEFLECTION  ANGLE  FIXED 

4  EFFECTOR  DEFLECTION  ANGLE  VARIABLE 

5  ACCIDENTAL  FAILURE  THRUST  LOSS 

6  ACCIDENTAL  FAILURE  EFFECTOR  REVERSAL 

7  PRINT 

8  HOLD 

9  RESTART  ADJUST  CONSTANTS  IN  PROGRAM 

10  RESTART  RETAIN  CONSTANTS 


Figure  17 
RTS5D  Modes  of  Operation 
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The  RESTART  mode  re-initialized  the  program  and  then 
provided  an  opportunity  to  the  operator  to  change  the 
magnitude  of  any  constant  in  the  program.   This  feature 
allows  the  operator  to  carry  out  program  modeling  adjust- 
ments without  the  delays  associated  with  re-compiling  a 
modified  digital  program,   Adjustments  to  constant  magni- 
tudes are  made  via  a  peripheral  teletypewriter  (TTY-35) 
(See  Fig.  5A) . 

The  ACCIDENTAL  FAILURE  FFFECTOR  REVERSAL  mode  deflected 
the  nozzle  angle  of  the  starboard  thruster  6.-  (See  Fig.  3) 
to  180**.   This  simulated  the  accidental  reversal  of  an 
outboard  thruster  at  full  power.   Simultaneously  the 
program  computed  and  displayed  the  performance  indices 
discussed  in  the  Software  section  of  this  report.   Exit 
from  this  mode  was  accomplished  when  the  pilot  had  maneuvered 
the  craft  back  on  track.   At  this  time  the  program  would 
"hold"  while  displaying  the  last  computed  values  of  all 
variables.   The  hold  condition  would  continue  until  a 
restart  was  accomplished  by  the  operator. 

The  ACCIDENTAL  FAILURE  THRUST  LOSS  mode  was  entered  by 
decreasing  the  thruster  setting  of  any  one,  or  combination 
of,  thrusters  (See  Fig.  5B) .   Exit  from  this  mode  was 
identical  to  the  accidental  failure  effector  reversal 
described  previously. 

The  FIXED  STEP  SIZE  INTEGRATION  mode  removes  the  real 
time  clock  input  and  utilizes  an  operator  specified  time 
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increment  on  each  iteration  of  the  program  loop.   This  is 
a  non  real  time  mode.   The  output  displays  and  other  modes 
are  functional  in  this  mode. 
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V.       RTS5D    RESPONSE    CHAR?^CTERISTICS 

A.  VALIDATION 

The  DBSIM5D  program  was  utilized  as  the  benchmark  to 
validate  the  RTS5D  simulation.   A  DBSIM5D  sequence  of  runs 
was  made  with  step  effector  deflection  angles  of  5",  10°, 
15**,  20°,  25°,  and  30°.   The  runs  were  initiated  with  the 
craft  at  56  knts  and  traveling  straight  ahead  just  prior 
to  effector  deflection.   The  first  peak  overshoot  value 
and  the  steady  state  value  attained  by  the  variables 
u  (surge)  ,  v  (sway)  ,  r  (yaw  rate)  ,  <^    (roll)  ,  and  6  (pitch) 
were  used  to  document  error  in  the  RTS5D.   The  results  of 
the  comparison  are  plotted  as  error  versus  effector  deflection 
angle  in  Fig.  18,  Fig.  19,  Fig.  20,  Fig.  21,  and  Fig.  22. 

B.  ACCIDENTAL  FAILURE  ANALYSIS 

The  response  of  the  RTS5D  when  subjected  to  accidental 
failure  mode  scenarios  was  examined.   Two  primary  situations 
involving  the  requirement  for  the  pilot  to  react  to  the 
failure  and  maneuver  the  craft  back  onto  its  original  track 
were  analyzed.   The  first  accidental  failure  situation 
consisted  of  reversing  the  starboard  thruster  T^q  (See 

Fig.  3),  such  that  6,^  equalled  180°  while  the  vehicle  was 
tracking  straight  ahead  at  56  knts.  The  second  situation 
consisted  of  again  reversing  the  T,q  thruster  to  6,^  equal 
180°  with  the  additional  requirement  that  the  failure  occur 
while  the  vehicle  was  established  in  a  steady  state  full  right 

turn  (6,=5o=(Sq=6t -.=-30°)  .  To  assume  that  the  "worst  case"  accidental 
/   o   9   10 
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Figure  18 
Effector  Angle  vs  Surge  Error 
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Figure  19 
Effector  Angle  vs  Sway  Error 
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Figure  20 
Effector  Angle  vs  Yaw  Rate  Error 
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Figure  21 
Effector  Angle  vs  Roll  Error 
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Figure  2  2 
Effector  Angle  vs  Pitch  Error 
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failure  scenario  prevailed  throughout  the  conduct  of  the 
two  failure  situations  described  above,  the  pilot  was  not 
allowed  to  reduce  the  thrust  magnitudes  of  either  the 
reversed  effector  T,q  or  the  remaining  thrusters  T_ ,  T- , 
and  Tg.   The  pilot  was  required  to  recognize  the  failure  had 
occurred,  respond  with  direct  real  time  man  generated  control 
efforts  via  the  "joystick"  (See  Fig.  5A)  and  maneuver  the 
vehicle  back  onto  the  "north  road". 

In  each  case  the  exercise  was  started  by  initializing 
the  craft  at  the  steady  state  conditions  for  either  a 
straight  ahead  run  or  a  -30**  effector  angle  right  turn. 
The  appropriate  steady  state  values  for  each  initialization 
used  are  shown  in  Table  IV. 


TABLE  IV.   STEADY  STATE  CRAFT  CONDITIONS  PRIOR 
TO  ACCIDENTAL  FAILURE  TEST 


VARIABLE 

STRAIGHT  AHEAD 

RIGHT 

TURN 

u 

87. 

ft/sec 

73.99 

ft/sec 

V 

0. 

ft/sec 

-7.58 

ft/sec 

R 

0. 

deg/sec 

2.9 

deg/sec 

<}> 

0. 

deg 

1.2 

deg 

e 

1. 

deg 

1.19 

deg 

^ 

0. 

deg 

0. 

deg 
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Termination  of  the  exercise  occurred  whenever  the  pilot 
simultaneously  satisfied  two  recovery  requirements: 

1)  The  craft's  navigation  coordinates  must  be  within 
50  ft  of  the  north  road  centerline  (+X    axis) . 

2)  The  translation  velocity  of  the  vehicle  in  the  Y„,„ 

NAV 

direction  must  be  less  than  plus  or  minus  1  ft/sec. 
It  was  considered  that  the  craft  was  completely  under  control 
once  the  vehicle  had  fulfilled  the  above  requirements.   It 
is  instructional  to  note  that  the  recovery  criteria  not 
only  required  that  the  pilot  counter  any  asymmetric  turning 
moments  and  sway  velocities  induced  by  the  failure,  but 
that  he  maneuver  the  vehicle  back  to  its  original  track. 
The  results  of  a  series  of  accidental  thrust  reversal  tests 
are  shown  in  Fig.  2  3A  and  Fig.  2  3B.   The  maximum  Y  ^^  (See 
Fig.  1)  translation  was  82.2  ft  for  the  straight  ahead  run 
case  and  571.2  ft  for  the  right  turn  effector  reversal  case. 
It  was  observed  that  in  both  cases  the  craft  retained  the 
capability  to  maneuver  back  to  the  original  track.   Further, 
the  corrective  action  required  in  both  cases  was  simply  to 
command  the  three  operational  effectors  to  the  full  left 
(+30")  position  and  hold  that  until  the  approach  to  the 
original  track,  at  which  time  the  effectors  were  modulated. 
Consequently  the  fastest  recovery  can  be  expected  to  come 
from  the  operator  who  is  quickest  to  recognize  that  an 
effector  reversal  has  occurred.   The  initial  optimal  control 
algorithm  in  both  cases  was  full  left  rudder. 


79 


100.0 


80.0 


60.0 


40.0  - 


20.0  - 


0.0 


5    6    7 
RUN  NUMBER 


10 


Figure  2  3A 
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Figure  23B 

[AV  Transla- 
Effector  Reversal,  Right  Turn  Run 
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Repeated  runs  made  with  the  same  operator  demonstrated 
a  very  steep  learning  curve.  Rapid  improvement  in  initial 
scores  was  demonstrated  in  every  case. 
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VI.   CONCLUSIONS 

The  ability  to  implement  a  real-time,  man-controlled 
simulator  exhibiting  the  steady  state  dynamics  of  the  3K-SES 
has  been  established.   By  requiring  the  right  hand  side  of 
the  simplified  equations  of  motion  to  be  restricted  to  a 
few  significant  terms,  that  portion  of  the  iteration  loop 
delay  dedicated  to  computation  and  the  AD/DA  process  has 
been  reduced  to  an  acceptable  magnitude  of  approximately 
35  msecs.   The  remaining  portion  of  the  total  iteration  time 
of  approximately  84  msec  was  required  for  the  subroutine's 
outputting  the  "alpha-numeric"  presentation  of  the  craft's 
states  onto  a  graphics  display  console.   This  low  value  of 
graphics  presentation  time  resulted  from  the  utilization 
of  time  multiplexing  of  the  data  to  be  displayed.   Using 
100  msecs  as  a  maximum  allowable  iteration  time  additional 
terms  for  the  right  hand  side  can  be  incorporated  to  generate 
a  more  accurate  simulation.   From  Section  III  on  rectuangular 
integration  error  analysis  it  was  concluded  that  the  RTS5D 
loop  delay  must  remain  less  than  100  msec  to  remain  credible. 
This  concurs  with  those  conclusions  reached  by  A.L.  Greensite 
in  Ref.  9.   Once  the  100  msec  loop  delay  point  has  been 
reached  further  expansion  of  the  right  hand  side  of  the  equa- 
tions of  motion  would  come  from  development  of  a  more 
efficient  graphics  system.   It  was  observed  that  while  the 
graphics  console  utilization  required  a  high  cost  in  terms 
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of  loop  delay,  the  interface  of  realistic  man-generated 
control  efforts  could  be  accomplished  in  no  other  available 
manner  due  to  limitations  on  existing  facilities. 

Results  of  operational  use  of  the  RTS5D  indicate  the 
craft  perspective  and  north  road  algorithm  developed  in 
Section  IV-D  and  used  to  generate  real  time  visual  cues 
for  the  pilot  to  be  extremely  effective  in  displaying  craft 
heading  changes  and  vehicle  displacements  in  Y  . 

The  results  in  Section  V-A  have  shown  the  RTS5D  to  be 
within  acceptable  error  limitations  of  less  than  2%  on  the 
average  when  compared  to  the  DBSIM5D  program  in  the  moderate 
speed  range  of  approximately  5  6  knts  and  over  the  full 
effector  deflection  range  of  plus  or  minus  30  degrees. 

The  results  in  Section  V-B  have  shown  that  the  3K-SES 
is  capable  of  being  controlled  satisfactorily  by  man  during 
accidental  failures  involving  both  loss  of  thrust  and 
effector  reversals.   From  these  results  it  was  shown  that 
the  optimal  control  effort  during  the  initial  stages  of  an 
effector  reversal  failure  was  simply  to  apply  full  opposing 
effector  angle  on  the  remaining  operational  effectors.   This 
input  was  held  by  the  pilot  until  the  craft  returned  to  the 
vicinity  of  the  original  track  at  which  time  the  operational 
effectors  deflection  angles  had  to  be  man-manipulated 
slightly  to  maintain  track.   Consequently  computer  assisted 
maneuvering  controls  would  be  somewhat  helpful  but  do  not 
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appear  to  be  a  hard  requirement  during  effector  reversal 
scenarios  in  the  3K-SES. 

It  was  also  concluded  that  the  RTS5D  program  can  be 
readily  expanded  to  a  full  6  degree  of  freedom  simulator 
with  the  addition  of  simplified  equations  representing  the 
heave  motion  characteristics  of  the  3K-SES. 
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VII.   RECOMMENDATIONS 

The  following  recominendations  are  made  for  future 
studies  utilizing  the  RTS5D  to  simulate  the  3K-SES: 

1.  The  current  right  hand  side  of  the  craft's  equations 
of  motion  require  expansion  to  include  more  terms  to 
accurately  represent  the  3K-SES  over  a  wider  operating  range. 

2.  Incorporate  the  sixth  degree  of  freedom,  heave, 
into  the  simulation. 

3.  Include  "sea-state"  input  whenever  more  computational 
time  becomes  available. 

4.  Program  the  propulsion  plant  dynamics  to  allow 
constant  power  as  an  input  in  lieu  of  constant  thrust. 

5.  Design  a  dedicated  micro-processor  or  mini-computer 
to  assist  the  operator  in  maneuvering  control  with  the  intent 
of  more  fully  examining  the  computer  assisted  controls 
requirement  in  the  3K-SES. 

6.  Incorporate  actual  3K-SES  control  hardware  into  the 
RTS5D  to  allow  more  realistic  "hands  on"  testing  of  not  only 
the  control  hardware  itself  but  also  the  operator  personnel. 
Currently  a  project  is  being  designed  to  test  operator 
response  under  various  command  and  failure  conditions. 
Results  will  be  reported  in  the  near  future. 
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APPENDIX  A 
RTS5D  PROGRAM  NOMENCLATURE 

All  A/D  trunk  500  voltage 

A21  A/D  trunk  501  voltage 

A22  Added  mass  coefficient  in  yaw  moment  equation 

A31  A/D  trunk  5  02  voltage 

A33  Added  mass  coefficient  in  roll  moment  equation 

A34  Added  mass  coefficient  in  pitch  moment  equation 

A41  A/D  trunk  503  voltage 

A51  A/D  trunk  504  voltage 

A61  A/D  trunk  5  05  voltage 

A71  A/D  trunk  5  06  voltage 

A81  A/D  trunk  507  voltage 

A91  A/D  trunk  510  voltage 

AASTOP  performance  index  limit 

ABSYD  1 Y  I 
'  o ' 

ABSYOD  I Y  I 
'  o  ' 

AD  array  of  A/D  lines 

AlCOA  speed  line  index  counter 

AM  mass 

APITCH  pitch  angle  in  degrees 

AP0TK5  voltage  on  pot  E  of  thruster  console 

APRINT  #  iterations  between  print  execute  commands 

AR  analog  turn  rate  voltage 
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ARATE 

turn  rate 

in  degrees 

ARDOT 

analog  r  voltage 

AROLL 

roll  angle 

k 

ASTOP 

performance  index  J,  limit 

AU 

analog  velocity  voltage 

AUDOT 

analog  u  voltage 

AV 

analog  sway  velocity  voltage 

AVDOT 

• 

analog  v  voltage 

AX 

X-Y 

plotter  scale  factor 

AXD 

X-Y 

plotte 

;r  X   coordinate 
o 

AXST 

initial  condition  X-coordinate  on  sea  track 

AY 

X-Y 

plotte 

ir  scale  factor 

AYD 

X-Y 

plotte 

;r  Y   coordinate 
o 

AYST 

initial  condition  Y-coordinate  on  sea  track 

BlOl 

D/A 

trunk 

T-432 

Bll 

D/A 

trunk 

T-420 

Bill 

D/A 

trunk 

T-433 

B21 

D/A 

trunk 

T-421 

B31 

D/A 

trunk 

T-422 

B41 

D/A 

trunk 

T-423 

B51 

D/A 

trunk 

T-424 

B55 

D/A 

trunk 

T-425 

B61 

D/A 

trunk 

T-426 

B71 

D/A 

trunk 

T-427 

B81 

D/A 

trunk 

T-430 

B91 

D/A 

trunk 

T-431 

BETA 

drift  angle 
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BOAT  subroutine  to  generate  perspective  of  SES 

CDP  bow  real  forces  lumped  coefficient 

CDX  drag  forces  lumped  coefficient 

CDY  sway  forces  lumped  coefficient 

CDZP  sidewall  rolling  moment  lumped  coefficient 

DAL  digital  to  analog  call 

DELT  iteration  loop  time 

DELTA  commanded  iteration  loop  time 

DTIME  V  time 

DTIMPLT  scaled  value  for  v  time  history 

DXPLOT  same  as  DTIMPLT 

DYPLOT  scaled  value  for  v 

DYREPET  restart  value  for  v 

EFLG  error  flag  for  TEXTO  call 

ETIME  r  time 

EYREPET  restart  value  for  r 

EXPLOT  scaled  r  time 

EYPLOT  scaled  r 

Fl  stern  buoyancy  force 

F2  bow  seal  pitch  force 

F3  bow  buoyancy  force 

FSP  port  sidewall  buoyancy  force 

FSS  starboard  sidewall  buoyancy  force 

FTIME  V  time 

FTIMPLT  scaled  v  time 

FXPLOT  same  as  FTIMPLT 

FYPLOT  scaled  v 
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FYREPET  restart  value  for  v 

G  gravitational  acceleration 

GRAPHO  subroutine  to  project  image  (non  alpha-numeric) 

GTIME  R  time 

GTIMPLT  scaled  R  time 

GXPLOT  same  as  GTIMPLT 

GYPLOT  scaled  R 

GYREPET  restart  R  value 

HEAD  craft  heading 

HOLD  subroutine  to  freeze  display /program 

HTIME  u  time 

HTIMPLT  scaled  u  time 

HXPLOT  same  as  HTIMPLT 

HYBETA  analog  B 

HYBRID  analog  flag  coefficient  (0  or  1) 

HYPLOT  analog  y  axis  sea  track  coordinate 

HYPSI  analog  heading  (radians) 

HYR  analog  turn  rate  (rad) 

HYRDOT  analog  r 

HYREPET  analog  restart  y  coordinate  on  sea  track 

HYSESX  analog  y  position  of  craft 

HYSESY  analog  x  position  of  craft 

HYU  analog  u 

HYUDOT  analog  u 

HYV  analog  v 

HYVDOT  analog  v 
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HYXO      analog  X 

• 

HYXODOT   analog  X 

^   o 

HYYO      analog  Y 

^   o 

* 

HYYODOT    analog  Y 

^   o 

ICO  multiplex  graphics  index 

ICOA  speed  line  index 

ICRAFT  vehicle  designator 

lER  error  flag 

IIPR  time  history  limit 

IIPRD  time  history  limit 

IIPRE  time  history  limit 

IIPRF  time  history  limit 

IIPRG  time  history  limit 

IIPRH  time  history  limit 

IIPRO  time  history  limit 

IIPRP  time  history  limit 

IIPRQ  time  history  limit 

IIPRR  time  history  limit 

IIPRS  time  history  limit 

IIX  print  counter 

IJ  craft  perspective  counter 

IJA  craft  perspective  counter 

IJAD  time  history  counter 

IJAE  time  history  counter 

IJAF  time  history  counter 

IJAG  time  history  counter 

UAH  time  history  counter 
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IJAO  time  history  counter 

IJAP  time  history  counter 

IJAQ  time  history  counter 

IJAR  time  history  counter 

IJAS  time  history  counter 

IJD  time  history  counter 

IJE  time  history  counter 

IJF  time  history  counter 

IJG  time  history  counter 

IJH  time  history  counter 

IJO  time  history  counter 

UP  time  history  counter 

IJQ  time  history  counter 

IJR  time  history  counter 

US  time  history  counter 

IL  graphics  console  display 

ILAl  graphics  console  display 

ILA2  graphics  console  display 

ILA3  graphics  console  display 

ILA4  graphics  console  display 

ILAS  graphics  console  display 

ILA6  graphics  console  display 

ILA7  graphics  console  display 

ILAS  graphics  console  display 

ILA9  graphics  console  display 

ILAIO  graphics  console  display 


array 
array 
array 
array 
array 
array 
array 
array 
array 
array 
array 
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ILAll  graphics  console  display  array 

ILA12  graphics  console  display  array 

ILA13  graphics  console  display  array 

ILA14  graphics  console  display  array 

ILA15  graphics  console  display  array 

IIiA16  graphics  console  display  array 

ILA17  graphics  console  display  array 

ILA18  graphics  console  display  array 

ILA19  graphics  console  display  array 

ILA20  graphics  console  display  array 

ILA21  graphics  console  display  array 

ILA22  graphics  console  display  array 

ILA2  3  graphics  console  display  array 

ILA24  graphics  console  display  array 

IIiA25  graphics  console  display  array 

ILA26  graphics  console  display  array 

ILA27  graphics  console  display  array 

ILA2  8  graphics  console  display  array 

ILA2  9  graphics  console  display  array 

ILAP  graphics  console  display  array 

IN  failure  mode  counter 

IPACK  subroutine  to  load  display  arrays 

IPERFT  performance  index 

IPLOT  graphics  console  display  array 

IPR  time  history  counter  limit 

IPRD  time  history  counter  limit 
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IPRE  time  history  counter  limit 

IPRF  time  history  counter  limit 

IPRG  time  history  counter  limit 

IPRH  time  history  counter  limit 

IPRO  time  history  counter  limit 

IPRP  time  history  counter  limit 

IPRQ  time  history  counter  limit 

IPRR  time  history  counter  limit 

IPRS  time  history  counter  limit 

ISWA  display  array 

ITEXT  subroutine  to  display  alpha-numeric  data 

IVEIWA  craft  perspective/north  road  array 

IVIEWAA  sea  track  array 

IVIEWB  sea  track  axis  array 

IVIEWC  u  time  history  array 

IVIEWD  V  time  history  array 

IVIEWE  r  time  history  array 

IVIEWF  V  time  history  array 

IVIEWG  r  time  history  array 

IVIEWH  u  time  history  array 

IVIEWO  p  time  history  array 

IVIEWP  q  time  history  array 

IVIEWQ  p  time  history  array 

IVIEWR  q  time  history  array 

IVIEWS  analog  sea  track  array 

IVIEWZ  time  history  axis  array 
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IX  X  axis  moment  of  inertia 

J  time  history  counter 

JD  time  history  counter 

JE  time  history  counter 

JF  time  history  counter 

JG  time  history  counter 

JH  time  history  counter 

JO  time  history  counter 

JP  time  history  counter 

JQ  time  history  counter 

JR  time  history  counter 

JS  time  history  counter 

L3  lever  arm  for  buoyancy  force 

LD  draft  of  sidewall 

LP  side  thrust  component  lever  arm 

LZ  lever  arm  for  sway  drag  force 

N  clock  interrupt  count 

NAD  AD  trunk  line  array 

NDA  DA  trunk  line  array 

00  lever  arm  for  thrust  side  component 

OTIME  p  time 

OTIMPLT  scaled  p  time 

OXPLOT  same  as  OTIMPLT 

OYPLOT  scaled  p 

OYREPET  restart  p  value 

P  roll  rate 

PB  plenum  pressure 
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PDOT  roll  acceleration 

PDOTM  max  roll  acceleration 

PHI  roll  angle 

PMAX  max  roll  rate 

POTK  thruster  console  array 

POTS  control  stick  array 

POTSA  function  box  array 

PSI  craft  heading 

PTIME  q  time 

PTIMPLT  scaled  q  time 

PXPLOT  same  as  PTIMPLT 

PYPLOT  scaled  q 

PYREPET  restart  q  value 

Q  pitch  rate 

QDOT  pitch  acceleration 

QDOTM  max  q 

QMAX  max  q 

QTIME  p  time 

QTIMPLT  scaled  p  time 

QXPLOT  same  as  QTIMPLT 

QYPLOT  scaled  p  value 

QYREPET  restart  p  value 

R  craft  turn  rate  r 

RDOT  r 

RDOTM  max  r 

READCLOCK  subroutine  to  sample  real  time  clock 
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RESET 

RHO 

RMAX 

RPERFl 

RPERF2 

RPERF3 

RTIME 

RTIMPLT 

RX 

RXPLOT 

RY 

RYPLOT 

RYREPET 

RZ 

SI 

S2 

S3 

S4 

SESX 

SESY 

SFl 

SF14 

SLIP 

SPDLIN 

SPEED 

STARTCLOCK 


subroutine  to  reset  analog  computer 

density  of  water 

maximum  yaw  rate  of  craft 

Jn 


q  time 
scaled  q  time 

same  as  RTIMPLT 

scaled  q  value 
restart  q  value 

I™ 


Tq  lever  arm 
T,Q  lever  arm 
Tg  lever  arm 
T_  lever  arm 


scaled  X 
o 

scaled  Y 
o 

scale  factor  in  time  history  arrays 
scale  factor  in  time  history  arrays 
Beta  influence  coefficient 
speedline  lower  screen  limit 
craft's  velocity  in  knts 
subroutine  to  start  clock 
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STIME  analog  time 

STIMPLT  scaled  analog  time 

SXPLOT  HYSESX 

SYPLOT  HYSESY 

SXREPET  restart  value  for  HYSESX 

SYREPET  restart  value  for  HYSESY 

T  total  thrust 

TIO  thruster  number  4 

TlOC  forward  vector  component  of  TIO 

TIOMAX  max  TIO  value 

TIOS  side  vector  component  of  TIO 

T7  thruster  number  1 

T7C  forward  vector  component  of  T7 

T7MAX  max  T7  value 

T7S  side  vector  component  of  T7 

T8  thruster  number  2 

T8C  forward  vector  components  of  T8 

T8MAX  max  value  of  T8 

T8S  side  vector  component  of  T8 

T9  thruster  number  3 

T9C  forward  vector  component  of  T9 

T9MAX  max  value  of  T9 

T9S  side  vector  component  of  T9 

TEXTO  subroutine  to  display  alpha-numeric  data 

TFORW  summation  of  forward  thrust  components 
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THETA  pitch  angle 

TIM  u  time 

TIME  time 

TIMPLT  scaled  TIM 

TINT  real  clock  interrupt  frequency 

TITLE  program  name  display  array 

TITLEO  program  name  display  array 

TITLEl  program  name  display  array 

TMAX  max  T 

TSIDE  sum  of  side  thrust  components 

TYAW  sum  of  yaw  moment 

U  forward  velocity  (surge) 

UDOT  u 

UDOTM  max  u 

UMAX  max  u 

UMUL  scale  factor  in  speedline  generation 

V  side  velocity  (sway) 

VCD  thruster  console  pot  array 

VDOT  V 

VDOTM  max  v 

VMAX  max  v 

VS  ■  total  velocity 

W  heave  velocity 

WDOT  heave  acceleration 

WE  width  of  bow  seal 

WRITEOXCK  subroutine  to  assign  clock  interrupts  value  to 
variable 
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WW        lever  arm  for  sway  drag  forces 

XO         X 

o 

• 

XODOT      X 
o 

XDX  boat  perspective  scale  factor 

XGX  boat  perspective  scale  factor 

XPLOT  display  array 

XREPET  reset  value  for  SESX 

XST  scale  factor  for  analog  plot 

YO         Y 
o 

• 

YODOT      Y 
O 

YHIGH  point  D  vertical  position  constant 

YPLOT  display  array 

YREPET  reset  value  for  SESY 

YST  scale  factor  for  analog  plot 

Z  operational  effector  angle 

ZIO  effector  number  four  angle 

Z2  upper  effector  angle  limitation 

Z3  lower  effector  angle  limitation 

Z4  zero  angle  upper  limit 

Z5  zero  angle  lower  limit 

Z7  effector  number  one  angle 

Z7I  modified  effector  number  one  angle 

Z8  effector  number  two  angle 

Z8I  modified  effector  number  two  angle 

Z9  effector  number  three  angle 
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Z9I  modified  effector  number  three  angle 

ZAB  scale  factor  for  ZIO 

ZI  modified  Z 

ZZZZ  fixed  step  effector  angle  input 
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APPENDIX  B 
RTS5D  Computer  Program  Listing 
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APPENDIX  C 
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